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Abstract 

We use the Gould-Hopper (GH) polynomials to investigate the Novikov- 
Veselov (NV) equation. The root dynamics of the a-flow in the NV equation 
is studied using the GH polynomials and then the Lax pair is found. In 
particulr, when N = 3,4, 5, one can get the Gold-fish model. The smooth 
rational solutions of the NV equation are also constructed via the extended 
Moutard transformation and the GH polynomials. The asymptotic behavior 
is discussed and then the smooth rational solution of the Liouville equation 
is obtained. 
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1 Introduction 



The Novikov-Veselov equation [3J H2J E7J E3] is denned by (U and t is real) 

U t = d 3 z U + B 3 U - 3d z (VU) -3d z (VU), (1) 

B z v = d z u. 

When z = z = x, we get the famous KdV equation (U = V = V) 

U t = 2U XXX - Y1UU X . 
The equation ([T]) can be represented as the form of Manakov's triad [23] 

H t =[H,A]-BH, 
where H is the two-dimensional Schrodinger operator 

H = d z B z + U 

and 

A = d 3 z -Vd z + B 3 z - Vd z , B = V Z + V- z . 
It is equivalent to the linear representation 

Htfj = 0, d t ip = Aijj. (2) 

We see that the Novikov-Veselov equation (OQ) preserves a class of the purely potential 
self-adjoint operators H. Here the pure potential means H has no external electric 
and magnetic fields. The periodic inverse spectral problem for the two-dimensional 
Schrodinger operator H was investigated in terms of the Riemann surfaces with some 
group of involutions and the corresponding Prym O-functions P EH ESI ESI l2"5] I5T]. 
On the other hand, it is known that the Novikov-Veselov hierarchy is a special 
reduction of the two-component BKP hierarchy [28l|33j(and references therein). In 
[28] , the authors showed that the Drinfeld-Sokolov hierarchy of D-type is a reduction 
of the two-component BKP hierarchy using two different types of pseudo-differential 
operators, which is different from the Shiota's point of view [31]. Finally, it is 
worthwhile to notice that the Novikov-Veselov equation (TTJ) is a special reduction of 
the Davey-Stewartson equation [20]. l2Tj . 

Let Hip = Hu = 0. Then via the Moutard transformation [U [231 12S] 

U (z, z) — > U (z, z) = U(z, z) + 2dd ln[z J (ipdco — udip)dz — {ipdu — uBip)dz] 
ip — > 6 = — / (ipduj — udip)dz — (ipdu — ujdip)dz, 
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one can construct a new Schrodinger operator H = d z d z + U and H\ = 0. The 
extended Moutard transformation was established such that [T?| 125] 

U{t,z,z) = U(t,z,z) + 2dd\n iW(ip,cu), 
V(t,z,z) = V(t,z,z) + 2dd\niW(ip,u) 

will also satisfy the Novikov-Veselov equation. Here the skew product W is defined 
by 



W(ip, oj) = j (ipduj — udip)dz — (ipdu) — udip)dz + [ipd 3 uj — ud 3 ip + ood 3 — tpd 3 uj 

+ 2(d 2 ipdu - dtpd 2 uj) - 2(d 2 i/)du - d^d 2 u) + 3V(ipduj - udip) 

- W(4)duj - udtl))}dt. (3) 

In particular, we can use U — V — as the seed solution. Then H = dd. Let us 
consider the holomorphic function P(z, t) and satisfy 

~dt ~ !h?' ( ' 

Then we have 
Theorem 1.1 f3^ 

Let Vi(t, z) and V2{t, z) be different holomorphic functions of z and satisfy (f^j. One 
defines u± — V\ + V\ and u 2 = V 2 + V 2 ■ Then 

U(t,z,z) = 2dd\niW(V 1 ,V 2 ), (5) 
V(t,z,z) = 2dd\niW(T 1 ,T 2 ) (6) 

is a solution of Novikov-Veselov equation. Here the skew product W is 

W(V 1 ,V 2 ) = V{P 2 - V 2 V X + J \{V[V 2 - V 1 V' 2 )dz + (V1P2 - VlV 2 )dz\ 



+ J [Pi'V 2 - V{P 2 + 2{V[V' 2 - V'[V' 2 ) + V{P 2 - Vi V 2 
+ 2{V 1 "V 2 - VtV 2 )]dt. (7) 



/// _ /// 



Moreover, ifVi(t,z) and V 2 (t,z) are polynomials, then the solution is rational in 
z, z, t. 

In [2], [Tl El H3] , the rational solutions and line solitons of the Novikov-Veselov 
equation ([I]) were constructed by the d-bar dressing method. To get these solutions, 
the scattering data has to be delta-type and the reality of U also puts some extra 
constraint on them. In these cases, the ^-function in ([5]) and can be expressed 
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determinant of some matrix. 
To study the dispersion relation 01]), the authors in [32J introduced the cr-flows 
for polynimial-type solutions 

F N {t, z) = z N + a x z N - x + a 2 z N - 2 + ■■■ + a N ^z + a N . (8) 

Then the flow (j3J) generates a linear flow 

a k = (N - K + 3)(N - k + 2)(N - k+l)a k ^ 3 , k = 1, 2, 3 • • • N. (9) 

It can be seen that 01,02 are conserved quantities. Indeed, 01, 2 , • • • , 0jv are the 
elementary symmetric polynomials in the roots qi, q 2 , ■ ■ ■ , qn of ¥ N (z): 

N 

= a 2(q) = ^qiqj, 

i=l i<j 

0-3(9) = ~ ^2 qiQjqk,- ■ ■ ,cr N (q) = (-l^q^- ■ -q N . (10) 

i<j<k 

The integrable (even linear) evolution of = (01, a 2 , ■ • ■ , 0jv) induces a dynamical 
system on the symmetric product S N C of the complex roots . We call such a dy- 
namical system on S N C a 0-system. 

From (J2J), we see that given two solutions Vi(t, z) and V 2 {t, z), one obtains a solu- 
tion of the Novikov-Veselov solution by a substitution of e lXl Vi(t, z) and e lX2 V2(t, z) 
into (j7|). The Ai and A2 defined here are real- valued constants. Therefore, to each 
pair of holomorphic solutions of PJ, we can get an (S 1 x S' 1 )-family of solutions to 
the Novikiv-Veselov equation [32J. 

The paper is organized as follows. In section 2, we describe the Gould-Hopper 
polynomials using the generating function and establish the recursive relation. In 
section 3, one studies the root dynamics of 0-flow and the Lax pair is constructed. 
Also, the asymptotic behavior is discussed. In section 4, the smooth rational solu- 
tions are found using the Gould-Hopper polynomials and their asymptotic behavior 
is investigated. Section 5 is devoted to the concluding remarks. 

2 Gould-Hopper Polynomials 

In this section, we introduce the Gould-Hopper polynomials and use them to get 
the solutions of @. To investigate the polynomial-type solutions of (jl]), inspired by 
the work in [11], one utilizes the Gould-Hopper polynomials [10]. The generating 
function of the Gould- Hopper polynomials P^{t,z) is 

00 \ j\r 
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Indeed, the Gould-Hopper polynomials P]y(t,z) has the operator representation 

J.2Q6 + 3f)9 +4Q12 

P N (t, z) = e«*z» = [1 + tdl + *A + IA + + . . .]^. 
Then one gets 

e *ef( e ^) = e ^+* 3 *. (11) 

We remark that in general the Gould- Hopper polynomials are defined by Pj^ (t, z) = 
e td ™ z N . Here we take m — 3. 

One notices that the Gould-Hopper polynomials P^(t,z) are characterized by 
(@| and Pjv(0, z) — z N . For example, 

P = 1, Pi = z, P 2 = z 2 , P 3 = z 3 + 6t, P 4 = ^ 4 + 24tz, 

P 5 = z 5 + 60i2 2 , P 6 = 2 6 + 120te 3 + 360t 2 

P 7 = z 7 + 210te 4 + 2520z, P 8 = z 8 + 336te 5 + 10080* V 

P 9 = 2 9 + 504tz 6 + 30240tV + 60480* 3 

p w = z 10 + 720tz 7 + 75600tV + 604800t 3 z (12) 



Actually, we have 

[ft/731 

£ fc Z 



W3] .fc.JV-3* 



= JflE il(w _3t)! 
= ATP^t,*) (13) 

From the operation calculus, one has 

(z + 3td 2 z )P N ^(t, z) = P N (t, z), N>1. 

Hence we yield the recursive relation 

P N (t, z) = zP N ^(t, z) + 3t(N -1)(N- 2)P N . 3 (t, z). (14) 

We can see that if we consider the equation (j3J) with the initial data of analytical 
function 

oo 

P(0,Z) = Y, a NZ N , 
N=0 

then the formal solution is 

oo oo 

P{t, z) = e td " a * zN = J2 a N p N{t, z). (15) 

N=0 N=0 
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The successive operation of the operator (z + Ztd\ ) on the solution (Tl5|) can help us 
construct more solutions of (jl]). For example, if P(0,z) = sinz, then we have, 

ja% ■ y _ ja'i sr^ V L ) „2Ar+i _ \^ I l ) p ( , v \ 
6 SmZ ~ 6 2^(2N + 1)\ Z -2^(2N + 1)\ P2N+1 ^ Z) 

= sm(z-t). 

The last equation uses (TTTT) . Hence 

(z + 3td 2 z ) N sin(z -t), N = 0, 1, 2, 3, 4, • • • 

are also solutions of (j4j). 
Remark: Let's define 

ip(X) = e Xz+xH - e ~ Xz - xH = 2sinh(\z + XH). 

Then y?(A) satisfies 



On expanding 

one has 
Actually, 



(p(X)zz = AV(A), <p(X)t = ¥>(A) 2 
00 

^(A) = 5>A 2l+1 , 



i=0 



1 > 0. 



r 2i+i n 

v 2i+l-3fc + fe 



, fc!(2i + l-3fc)! 

It is known that <pi can be used to construct the Wronskian solutions of the KdV 
equation. The details can be found in [22J. 



3 Root Dynamics of cr-flows 



In this section, one uses the Gould-Hopper polynomials to study the root dynamics 
of the (x-flows ©• 
Let's write flSJ) as 

F N (t,z) = (z- qi (t))(z-q 2 (t))---(z-q N (t)). 

Then from the equation (j3J), one gets the root dynamics 

N 1 

<ij = ~ 6 5^ 7 w r ( 16 ) 

, (Qj-Qm )(qj-q n • 

m<n, jf=m,n J J 



For example, when N=3, we have 

g'i = -6 
q 2 = -6 
q 3 = -6 

For N=4, we get 

?i = " 6 b ^ v + 
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1 


- 9s) 


(93 


- 9i)(9s 
1 


- 92) 


(9i 


- 93) (9i 
1 


- 94) 


(92 


- 93) (92 
1 


- 94) 


(93 


- 9i) (93 
1 


- 94) 


(94 


- 9i)(94 


- 92) 



42 - " 6[ (92 - 9i)(92 - 93) + (92 - 93)(92 - 94) + (92 - 9i)(92 - ^ 

* " -><^^ + (^« + (^^ 

We notice that since o\ and a 2 are conserved quantities, one knows that 

N N 

1=1 8=1 

are conserved densities of (|T6|) . 

Now, we can investigate the properties of the root dynamics fTT6l) by the Gould- 
Hopper polynomials: 

• Initial Value Problem : The root dynamics of ex-flow can be solved by 

F N (t,z) = (z- qi (t))(z-q 2 (t))---(z-q N (t)) 

= P N (t,z) + C 1 P N _ 1 (t,z)+--- + C N P (t,z), (17) 

where the constants Ci, C 2 , • ■ ■ , Cjv-i, Cat are determined by the initial values 
of qi(0), g 2 (0), • • • , 9iv(0), that is, 

AT 

Ci = -X)ft(0), C 2 = 5^(6(0)^(0), 

i=l i<j 
i<j<fc 

= (-l)%(0)g 2 (0)---g w (0). 

Therefore, it is seen that the solutions qi(t), g 2 (t), - • • , (/jv(t) can be obtained 
algebraically 
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Lax pair: 

Firstly, we study the root dynamics of the Gould-Hopper polynomials, which 
correspond to the initial values (?i(0) = q 2 (0) = • • • = <?Ar(0) = 0. 
Let's define the N x N matrix by 



{0^+1 = 1, if i = 1,2,3, ~-,N ; 

a M _ 2 = -3t(i- l)(i-2), if i = 3, 4- • • , N - 1; 
0, otherwise. 

Then from the recursive relation (|T4j) . one knows that 

Pjv(t,z) = det(zI N - X(t)). 

For example, when iV = 3, 

X(t) = 



N=4, 




X{t) 



N=5, 



X(f) 





/ 





1 



















1 











-6t 








1 


) 




V 





-1* 


it 


o) 




/ 





1 










o \ 












1 










St 










1 










-1? 


it 








1 


V 










-36t 





o / 



(18) 



We can write X(t) as 

X(t) = R(t)QR-\t), 
where Q = diag(qi(t), q 2 (t), • ■ ■ , ?jv(t) and 



R(t) 



( P (qi,t) P (q 2 ,t) P (q 3 ,t) ••• P (q N ,t) \ 

Pi(qi,t) Pi(q 2 ,t) Pi(q 3 ,t) ■■■ Pi{q N ,t) 

P 2 {q u t) P 2 {q 2 ,t) P 2 {q 3 ,t) ••• P 2 (q N ,t) 

\P N (q u t) P N (q 2 ,t) P N (q 3 ,t) ■■■ P N (q N ,t)J 



(19) 
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For instance, when N 



and N=4, 



N=5, 



R(t) 



R(t) 



1 1 1 

Qi Qi Q3 

qI ol qI 



I 


1 


1 


1 


1 \ 




Ql 




Q3 


Qi 




Ql 


Ql 


qI 


Ql 


\ 


qf + Qt 


ql + Qt 


ql + 6t 


ql + Qt ) 



R(t) 



( 1 


1 


1 


1 


1 \ 


Qi 


Q2 


Qa 


Qa 


Qb 


Ql 


Ql 


Ql 


Ql 


Ql 


ql + Qt 


q\ + Qt 


4 + Qt 


ql + Qt 


qI + Qt 


\q\ + 24t 


qi + 24t 


ql + 24t 


qi + 24t 


Ql + 24t / 



From the initial value problem ( II 7p . we notice that the polynomials t n can be 
replaced by the elementary symmetric polynomials of the roots qi, q2, ■ ■ ■ , qN- 
Hence one has R(q). It can be seen that 



X(t) = RLR- 



where 



L = Q + [M,Q], M = R~ 1 R. 
For example, when N=3, 

/ 

L(t) = 



Qi 

^ x 93-92 
\ ^ 1 92-93 



Qi 



g2-93 

93-91 
Q2 

^91-93 



Q3 



Q3 



93-92 

92- 91 

93- 91 
91-92 

Q3 



I 



and N=4, using the Maple software, 



/ 



L(t) 



9i (91 



Qi 

-93)(9i-94)+6 



92(92-93)(92-94)+6 



(92-93)(92-94) 
9l(9l-92)(9l-94)+6 

(92-93)(94-93) 
9l(9l-92)(9l-93)+6 

(94-93)(94-92) 



(91 



g'2 (g2 



-93)(9i" 

Ql 

-94) (92- 



-94) 



9i)+6 



Since 



X(t) 



dX(t) 
dt 



0. 



(94 -93)(9l-93) 
92(g2-93)(92-9l)+6 

(94 -93)(94 -9l) 



-3(i-l)(i-2), 



93(93-g2)(g3-g4)+6 

(gi-g2)(gi-g4) 
93(g3-gi)(g3-g4)+6 

(92-g4)(g2-gi) 
Q3 

93(g3-gi)(g3-g2)+6 
(94-gi)(g4-g2) 



9 4 (94-92)(g4-g3)+6 \ 

(gi-92)(gi-g3) 
94(94-gi)(g4-g3)+6 

(92-93)(g2-gi) 
94(94-gi)(g4-g2)+6 



(92 



-93) (91- 
Q\ 



93) 



if i = 3,4- --,N- 1; 
otherwise, 
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we know X(t) is a nilpotent matrix and hence L is a nilpotent one, too. So 
tr(L r ) = tr[X(t)} r = 0, r = 1, 2, 3, . . . , .. 

Actually, a simple calculation yields 

L [f]+i = 0; Ar > 3; 

where [a] means the integer part of a. Now, 

d 2 X(t) 
dt 2 

will imply the Lax equation 

dL(t) 







dt 



[L, M). (20) 



For N = 3,4,5, we see that, by the Maple software, qi satisfies the following 
Goldfish model [5], a limiting case of the Ruijesenaars-Schneider system: 

& = 2 £_ft?L_. (21) 

The reason is that Pj, z = 3, 4, 5 are linear in t- variable (see the appendix). For 
iV = 6, we have from the diagonal terms of the Lax equation f[2"0"j) 

QiQi 5^7=1 ( some quadratic terms of q)qj + 720 
Qi — 2 > 1 s . 

Secondly, we consider the general case. Let's define the 2D Appell polyno- 
mials G n (z,t) by means of the generating function [4]: 

A™ 



A{\)e Xz+xH = VG„( 2i ^, (22) 
' n! 

n=0 

where 

r^.s being constants and Tq — 1. Then one has the following formula, noting 
thatffl = (/_J, 

- / N \ 



ITJV — 

h=0 



N-h 
10 



It's easy to see that the polynomials G n (z,t) also satisfy the linear equation 
OH). When comparing ( 1T7|) with (|2"51) . we have 

T N ~h = / AT \ • ( 24 ) 
\N-h) 

Now, it's suitable to introduce the coefficients of the Taylor expansion 

A'(\) _ ^ A™ 

K 1 n=0 

It can be seen that the coefficients ot n can be expressed by r , Ti, • • • , r n+ i 
(or the initial values fl24"]) ). For examples, 



a = r 1; ai = r 2 -r?, a 2 = r 3 - 3r x r 2 + 2Tl 
« 3 = r 4 + i2r?r 2 - 4r x r 3 - 3r| - 6T\, 



The recurrence relation for the 2D Appell polynomial Gjv(2, t) can be written 
as follows [1]: 

G (z,t) = 1 

G N {z,t) = {z + a )G N _ 1 {z,t)+3t{N -1){N -2)G N _ 3 {z,t) 
N ~ 2 /N-l\ 

+ J2[ k J^N-k-iG k (z,t). (25) 

A simple calculation can yield 

G (z,t) = 1, Gi(z,t) = z + a , G 2 (z, t) = (z + a ) 2 + a\ 

G 3 (z,t) = (z + a ) 3 + 3ati(z + a ) + a 2 + Qt 

G 4 (z, t) = (z + a ) 4 + Ga^z + a ) 2 + (4a 2 + 2U){z + a ) + «s + 3a 2 

G 5 (z,t) = (z + ato) 5 + I0a 1 (z + a ) 3 + (I0a 2 + 60t)(z + a ) 2 

+ (5a 3 + 15a?) (z + a ) + 60tai + 10aia 2 + a 4 . 

When A(X) = 1, this recursive relation becomes f[T41) . Hence the relation f[2"5]) 
is a generalization of (j!4p for arbitrary initial data. The matrix corresponding 
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to (TTSl) can be constructed as follows: 



X(t) 



a i,i+l 


= 1, 






0-1,1-1 




0>i,i-2 




/ Q>i,i— 3 




&i,i—k 





J)ai, 



)< 

) (6t + 

l) a 3, 

(i-fc-lH> 



3 
i-4 



if i 


= 1,2,3,- 


••,JV-1 


if i 


= 1,2,3,- 


••,JV; 


if i 


= 2,3,4,- 


•-,iV; 


if i 


= 3,4,5,- 


■-,iV; 


if Z 


= 4,5,6,- 


-.,N; 



ifi = k + l,k + 2,k + 3, 



,N; 



Similarly, one has 

GatOM) = det(zI N - X(t)). 
For instance, when A^ = 5, we get 

/ -a 1 

— «i — «o 1 

AT(t) = -(6t + a 2 ) -2«i -a 

— 0:3 — (18i + 3«2) — 3ai 

\ — «4 — 4«3 — (36t + 6«2; 

Also, one can write AT(t) as 





1 

-Qt 

-4«i 






1 

-a J 



X(t) = R(t)QR-\t), 

where Q = diag(qi(t), q^it), ■ • ■ , <Zjv(t) and R(t) is defined as ffT9l) with P m (qi, t) 
being replaced by G m (g«,t). Then one follows the previous procedures and 
finally can get the Lax equation (|2"U1) for general case. Therefore the root 
dynamics (TIB"]) is Lax-integrable. But the computations are more involved and 
one doesn't pursuit them here. 

We notice here that for A^ = 3,4,5 the root dynamics of G^ also satisfies 
the Gold-fish model (EH). 



Asymptotic behavior 

It is known that the Gould-Hopper polynomial P^(t.z) has the scaling prop- 
erty: 



P N {t,z) 



N ~ , Z 
t~Pj 



t X/3- 



(26) 



where Pn{v) is the so-called Appell polynomials [I] (and references therein) 
in rj = ~t/3. For example, 

-21 



Pa(t, z) 



z 8 + 336te 5 + 1008(WV = tsfo 8 + 336r? 5 



10080rf 



tzPsiv)- 
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Then the k-th zero Xffl of Pn(v) determines the dynamics of the root q k , i.e., 

q k (t) = 

Since Pn(^X^) = 0,£ 3 = 1, one knows that the roots q k are located on the 
circles in the plane with time dependent radius or fixed at the origin. Finally, 
from the Initial value Problem ( TlT)) and fl26|) . we know that when t — » oo and 
z — > oo such that \z\ 3 /t — > constant, Pjv(i, z) plays the dominant role. Then 
one yields 

q k (t) -> 

Consequently, the roots asymptotically will follow diagonal lines. 



4 Smooth Rational Solutions of Novikov-Vaselov 
equation 

In [32], the blow-up solution of the Novikov-Vaselov equation via the extended 
Moutard transformation (J2J) is constructed. In this section, we establish smooth 
rational solutions for all time by the Gould -Hopper polynomials (|12|) . Some calcu- 
lations below need the Maple software. 



Example 1 

Let V\ = z 2 + z + 1 and V2 = —iz 2 — 2iz. Then a simple calculation can yield 
the imaginary part of W in (J7|) 

g 

M(x, y, t) = (x 2 + y 2 ) 2 + -a; 3 + Axy 2 + Ax 2 + Ax + At + 100 (27) 

From f l27|) . we can see that 

M(x, y, 0) (x 2 + y 2 ) 2 near r = \J ' x 2 + y 2 = 00. 

It can be verified that M(x, y, 0) is positive for all M. 2 . Also, M(x,y,t) is 
positive for all M 2 at any fixed time t > . Then the solution {7 in §5§ of the 
Novikov-Veselov equation is 

M 2 ' 

where 

Mi = -12(294 + 600a; 2 + 588?/ 2 + 8a; 3 + 888x + 12t + 3x 4 -6a; 2 ?/ 2 -2a;V 
+ 24x 2 t-3?/ 4 x + 24?/ 2 t + 36a;t-3?/ 4 + a; 5 ) 
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and 



M 2 = (3 x 4 + 6 x 2 y 2 + 3 y 4 + 8 x 3 + 12 xy 2 + 12 x 2 



+ 12x + 12t + 300) 



2 



At fixed time t , one knows U decays like 4? for r — > oo. Also, C/ tends 
asymptotically to zero at the rate t at any fixed point (x, y) when t — > oo. 

Example 2 

Let V\ = (z 3 + 6t) + 2iz and P 2 = —i(z 3 + 6t) + z. Then a simple calculation 
can yield the imaginary part of W in (j7j) 

f(x,y,t) = {x 2 +y 2 ) 3 +Ax 3 y + 8xy 3 + 2(x 2 +y 2 ) + 6t(2x 3 -Qxy 2 -y) + 3Qt 2 + 6000 

(28) 

From (l28~D. we can see that 



f(x, y, 0) ~ (x 2 + y 2 ) 3 near r = \/ x 2 + y 2 



oo. 



It can be verified that f(x,y,0) is positive for all M 2 . Also, letting ([28]) be 
equal to zero, one has 

t = (l/2)xy 2 + (l/12)y-(l/Q)x 3 

± (1/12) v^xV - 20xy 3 - 36xV - V - 20x 3 y - % 6 - 8x 2 - 23996, 

A simple calculation shows that the equation inside the square root is negative 
for all R 2 . Hence f(x,y,t) is positive for all R 2 at any fixed time t > 0. Then 
the solution U in (jSJ) of the Novikov-Veselov equation ([T|) is 



where 



Fx = ~[2Ax 7 y + 16x 6 + 24xV+ (432 + 216000 + 48 y 2 ) x 4 

+ (-24 y 5 - 144 1 + 48 y) x 3 + (432000 y 2 + 864 1?/ 3 - 96 y 4 ) x 2 
+ (-24 y 7 - 48 y 3 + 432 ty 2 + (432000 + 1728 1 2 ) y) x + 48000 
+ 432 yh - 32 y 6 + 216000 y 4 + 252 1 2 ] 

and 

F 2 = (x 6 + 3xV + 3xV + y 6 + 4x 3 ?/ + 8x?/ 3 + 2x 2 + 2?/ 2 
+ 12 tx 3 - 36 txy 2 - 6 ty + 36 1 2 + 6000) 2 

At fixed time t , one knows U decays like \ for r — >■ oo. Also, U tends 
asymptotically to zero at the rate at any fixed point (x,y) when t — > oo. 
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Example 3 

Let V\ = (z 4 + 24tz) + 2iz and V% = —i(z A + 24tz) + z. A tedious calculation 
shows that the imaginary part of W in ([7]) is 



y, t) = (x 2 +y 2 f+Qx 4 y+12x 2 y 3 -(18/5)y b +2x 2 +2y 2 +2At(x 2 +y 2 )(2x 3 -Qxy 2 +2At)+100 

Similarly, g(x, y, t) is positive for all M 2 at fixed time t > 0. Then the solution 
{7 in (0) can be written as 

Gi 



where 



.r 2 



G x {x,y,t) 

= -- [240 y 11 - 1680 xV - 1980 + (34560 tac - 5280 x 4 ) 

+ (160000 - 2160 x 2 ) y 6 + (480 + 218880 tx 3 - 3360 x 6 + 138240 t 2 ) 

+ (480000 x 2 - 9000 x 4 ) y 4 + (-120000 + 1200 x 8 + 57600 tx 5 

+ (-1382400 1 2 - 4800) x 2 ) y 3 + (3600 x 6 - 5760000 tx + 480000 x 4 ) y 2 

+ (1200 x 10 + 57600 x 7 t + (2400 + 691200 t 2 ) x 4 + 360000 x 2 ) y 

+ 20000 + 160000 x 6 + 1920000 tx 3 + 5760000 1 2 + 900 x 8 ] 

and 

G 2 (x,y,t) 

= (5 x 8 + 20 xV + 30 xV + 20 xV + 5 y 8 + 30 x 4 y + 60 x 2 y 3 - 18 y 5 + 10 
+ 10 y 2 + 240 tx 5 - 480 txV + 2880 x 2 t 2 - 720 ty A x + 2880 y 2 t 2 + 500) 2 

In this case, at fixed time, one knows that U decays like \ for r — > 00 ; 
moreover, [/ also tends asymptotically to zero at the rate 4 at fixed point 
(x, y) 7^ (0, 0) when t — > 00. We notice here that at (0, 0), U approaches 00 as 
t — )■ 00. 

Example 4 

Let Vi = 5 + 60te 2 ) + 2iz and V 2 = -i(z 5 + QOtz 2 ) + z. Then the imaginary 
part of in ([ID is 

h{x,y,t) = (^(x 2 + y 2 ) 5 -^x 5 y+^x 3 y 3 +(x 2 + y 2 )(2 + l2x 3 y-l2xy 3 ) 

+ 120 1 (x 2 + y 2 ) (x 5 - 3 y A x - 2 xV ) + 3600 1 2 (x 2 + y 2 f 
+ 20t(ll?/ 3 + 3x 2 ?/) + 120t 2 + 1000) 
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Similarly, h(x, y, t) is positive for all M 2 at fixed time t > 0. Then the solution 
U in P is 

Hi 
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where 



Hi(x,y,t) 

= -^[360xy 13 + (360 x 3 + 14400 1) y 11 - 144 y 10 + (14400 x 2 t - 1800 x 5 ) y 9 

+ (-720 x 2 + 900000 + 108000 1 2 ) y 8 + (237600 tx 4 - 777600 t 2 x - 3600 x 7 ) y 7 
+ (3840 x 4 + (3600000 + 432000 t 2 ) x 2 + 83520 tx) y 6 + (165600 tx 6 + 720 x 

- 10368000 1 3 - 1800 x 9 - 345600 t 2 x 3 ) + (-1120 x 6 + (5400000 

+ 648000 1 2 ) x 4 + 12000 tx 3 + (-129600000 1 - 15552000 t 3 ) x - 810000 1 2 ) y 4 
+ (360 x 11 - 50400 xH + 950400 x 5 t 2 - 2160 x 3 + 20736000 t 3 x 2 + (-216000 1 2 

- 1800000) x - 5760 1) y 3 + (2380 x 8 + (3600000 + 432000 1 2 ) x 6 - 25200 tx 5 

+ (-10368000 1 3 - 86400000 1) x 3 + 453600 x 2 t 2 + 62208000 1 4 + 518400000 1 2 ) y 2 

+ (360 x 13 + 21600 x 10 t + 518400 x 7 t 2 + 720 x 5 + 31104000 t 3 x 4 + (1920000 

+ 230400 1 2 ) x 3 + 17280 x 2 t + 1555200 1 3 + 12960000 1) y + 476 x 10 

+ (900000 + 108000 1 2 ) x 8 + 25200 x 7 t + (43200000 1 + 5184000 1 3 ) x 5 

+ 226800 xH 2 + (62208000 1 4 + 518400000 1 2 ) x 2 + 72000 + 8640 1 2 )} 

and 

H 2 (x,y,t) 

= (3 x 10 + 1 5 xV + 30 x 6 y A + 30 y 6 x A + 15 y s x 2 + 3y 10 + 26 x 5 y + 20 x 3 y 3 

+ 6 x 2 + 6 y 2 - 36 xy 5 + 360 x 7 t - 1800 tx 3 y 4 - 360 tx 5 y 2 - 1080 txy 6 

+ 10800 x 4 t 2 + 21600 t 2 x 2 y 2 + 10800 t 2 y A + 660 ty 3 + 180 tx 2 y + 360 1 2 + 3000) 2 

In this case, at fixed time, U decays like \ for r — > oo ; moreover, we see that 

-240(x 2 + y 2 ) -2A0zz 

U — > 7 ', — 7, TTT = 7 777 OS H OO I 29 I 

[30(x 2 + y 2 ) 2 + l] 2 [30z 2 z 2 + l] 2 V ; 

at fixed point (x, y), which is a stationary solution of (JTJ for 

3600 5V - 120 z 2 



2 



V 

{30z 2 z 2 + 1) 
We notice that if we define 

u(z, z)—lnU and V(z, z) = t^T' 
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then the stationary equation of the Novikov-Veselov equation ([TJ) will become the 
Tzitzeica equation JHJ [TH] 

Uz - = e u + ee~ 2u , (30) 

where e is an arbitrary constant. It can be verified that (|29i) satisfies the e = case, 
i.e., the Liouville equation, whose real solutions are given by 



-2k 



dS\2 



ln (TT#p (31) 

where S(z) is a locally univalent meromorphic function in some domain and k is 
a constant. For the solution corresponding to f )29|) . one knows that S(z) = z 2 and 
k = 30. 

For general case, let's choose 

V\ = Pjv(^ z ) + 2zz and P 2 = — iPN(t, z ) + 

We can expect the imaginary part of W in (j7J) is positive for all K 2 at any fixed 
time if we choose an appropriate constant. And the solution U(x,y,t) at any fixed 
time decays like piTT for N > 3. For N = 1, one refers to [13] ■ It could be 
interesting to study the asymptotic solutions when t — > oo and then one obtains the 
solutions ( 13~T|) of the Liouville equation or the ones of the Tzitzeica equation (l30]h 
One remarks here that in these examples, for each potential U(x, y, t) there exist 
infinitely many wave functions, which can be constructed by the Pfaffian [TJ of 
linear combinations of the Gould- Hopper polynomials. Hence it could be worthwhile 
to investigate the corresponding wave functions when t — > oo. But the computations 
are more involved and the details will be published elsewhere. 



5 Concluding Remarks 

In this paper we have studied the Novikov-Veselov equation using the Gould-Hopper 
polynomials. Firstly, one investigates the root dynamics of the so-called <r-flows and 
gets the Lax pair; moreover, one finds that when iV = 3, 4, 5, the root dynamics 
satisfies the Gold-Fish model. Although the Lax pair is established, only two con- 
served densities are found. The reason is that the Lax operator is nilpotent. Also, 
the asymptotic behavior is studied. Secondly, we construct smooth rational solu- 
tions using the Gould-Hopper polynomials and the skew product ([7]); besides, the 
asymptotic behaviors of the smooth rational solutions are discussed. 
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6 Appendix 



A Gold-Fish Model 

For the Goldfish Model 



QiQj 



■ Qi-Qj 



its initial value problem can be solved by the statement: z = q%{t),i = 1, 2, • • • , N 
are the N roots of the equation [5] 

gi(0) 1 

< y. — n l 



i=l 



It can be seen that it is a polynomial in z with coefficients linear in t. Then the 
special choices of initial datum can get the solutions of the root dynamics ffT6l) for 
the cases N = 3,4, 5. To illustrate it, we take N = 3 as an example. When N = 3, 
we have 

gi(Q) g 2 (Q) g 3 (Q) _ i 

z- qi (0) z-q 2 {0) z-q 3 (0) t 
After some calculations, one yields 

z 3 - z 2 [( gi (0) + ^(0) + g 3 (0)] + z[ qi (0)q 2 (0) + g 2 (0)g 3 (0) + gi(0)g 3 (0)] - gi(0)g 2 (0)g 3 (0) 
= tz 2 [ qi (0) + g 2 (0) + q 3 (0)\ - te[ ?1 (0)(g 2 (0) + g 3 (0)) + g 2 (0)( gi (0) + g 3 (0)) 
+g 3 (0)(gi(0) + g 2 (0))] + t[gi(0)g 2 (0)g 3 (0) + g 2 (0)gi(0)g 3 (0) + g 3 (0)gi(0)g 2 (0)]. (A.l) 

On the other hand, from (TlTj) . one knows qi(t), g 2 (t), qs(t) are the roots of the poly- 
nomial 

F 3 (z, t) = z 3 + 6t + C lZ 2 + C 2 z + C 3 

or 

z 3 + C lZ 2 + C 2 z + C 3 = -6t. (A.2) 
Comparing flA.lj) with flA.21) . we are able to get the following linear equations for 

gi(o),g 2 (o),g 3 (o) : 

gi(0) + g 2 (0) + g 3 (0) =0 

gi(0)(g 2 (0) + g 3 (0)) + g 2 (0)(gi(0) + g 3 (0)) + g 3 (0)( gi (0) + g 2 (0)) = 
gi(0)g 2 (0)g 3 (0) + g 2 (0)g!(0)g 3 (0) + g 3 (0)g x (0)g 2 (0) = -6. (A.3) 

So if the determinant of the matrix 

1 1 1 

g 2 (0) + g 3 (0) gi(0) + g 3 (0) g x (0)+g 2 (0) | (A.4) 
g 2 (0)g 3 (0) gi(0)g 3 (0) gi(0)g 2 (0) 
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is not equal to zero, then the initial velocities gi(0), ^(O), ^3(0) can be uniquely 
expressed by the initial positions ?i(0), g 2 (0), (^(O). For the cases N = 4,5, the 
linear equations (IA.3[) and the matrix (IA.4|) can be obtained similarly. 
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